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USING DYNAMIC INTERFEROMETRIC SYNTHETIC APERATURE RADAR 
(INSAR) TO IMAGE FAST-MOVING SURFACE WAVES 

[0001] The United States Government has rights in this invention pursuant to 
Contract No. W-7405-ENG-48 between the United States Department of Energy and 
the University of California for the operation of Lawrence Livermore National 
Laboratory. 

RELATED APPLICATION 
[0002] This application claims the benefit of U.S. Provisional Application No. 
60/421,423, filed October 24, 2002, and entitled, " Using Interferometric Synthetic 
Aperature Radar (InSAR) To Image Seismic Surface Waves," which is incorporated 
herein by this reference. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

[0003] The present invention relates to synthetic aperture radar (SAR). More 
particularly, the invention relates to a sampling method and apparatus for imaging 
fast-moving surface displacements in real time, such as seismic waves, utilizing both 
single-pass and repeat-pass synthetic aperture radar interferometry (InSAR). 
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State of Technology 

[0004] Seismology plays a dominant role in monitoring global natural (e.g., 
earthquakes) and man-made activity (e.g., underground mining). Seismology, 
however, is an inherently blind methodology, whereby the seismic waves are 
recorded by the effect they have on instruments as the seismic waves pass by. In 
addition, the type and quality of information available from seismology has been 
and continues to be limited by the relatively sparse spatial sampling on the earth's 
land surface and nearly complete lack of spatial sampling over the oceans. If these 
dynamic waves can be imaged in real time, actual movies of these waves can be 
produced yielding more information about their source and the structure they 
propagate through. 

[0005] InSAR (Interferometric Synthetic Aperture Radar) is a new geodetic 
standard providing spatially continuous maps of surface displacements to the 
millimeter level, with 25-meter horizontal (pixel) resolution. Operational benefits of 
InSAR include all weather, day or night, spatially continuous, wide area (100 by 100 
km image frames) and global coverage from space. While InSAR has predominantly 
been used to map and measure static deformations, the possibility exists for 
applying InSAR techniques to image dynamic deformations, such as those from 
seismic waves propagating along the earth's surface. 

[0006] Three typical ways of acquiring InSAR data are: 1) cross-track single pass 
(dual anterma), 2) cross-track repeat pass (single anterma), and 3) along-track single 
pass (dual antenna) . The single-pass methods require two SAR antermas moimted 
on the same platform for simultaneous acquisition. The repeat-pass method requires 
only one antenna that acquires data over the same area twice via a repeat pass, but 
with a slightly different viewing geometry. This approach requires precise 
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knowledge and predictability of the flight path and hence it is best suited to space- 
borne systems. 

[0007] Along-track interferometry (ATI), which is applicable to the remote 
measurement of velocities of objects moving on the groimd, images a ground region 
twice at slightly different times by producing two SAR images. The movement of the 
measured object in the images, e.g., water currents, causes a phase shift between the 
corresponding pixels in the SAR images. The moving surface produces a Doppler 
shift relative to the other stationary ground pixels,. A Doppler image can thus be 
created whereby only moving objects are visible in the image. Measurements from 
an inertial navigation system may be used to correct for unwanted phase shifts 
caused by aircraft yaw, pitch and roll. While a SAR is usually operated as a single- 
channel SAR, its antenna may be divided in the direction of movement and operated 
as a two-channel ATI SAR (i.e., an along-track interferometer). 
[0008] Background information on a single SAR line utilizing a single pass is 
disclosed in U.S. Patent No. 5, 945, 937, titled "ALONG-TRACK 
INTERFEROMETRIC SYNTHETIC APERTURE RADAR," issued August 31, 1999, 
to Fujimura, including the following: "An along-track interferometric SAR of the 
present invention includes a single SAR channel for acquiring SAR data by 
observing a target only once. A reconstructing section reconstructs two SAR images 
deviated in time from the SAR data to thereby output two reconstructed SAR 
images. A detecting section calculates a phase difference between the two 
reconstructed SAR images to thereby detect the velocity of the target in the SAR eye 
direction. The reconstructing section bisects, based on the fact that the frequency of a 
wave returned from the ground undergoes a Doppler shift due to the movement of 
a SAR relative to the groimd, the SAR data in the frequency domain and processes 
the bisected SAR data." 
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[00091 Despite existing conventional systems and methods, a need exists for an 
improved dynamic InSAR for real time assessment of surface displacements, such 
as, but not limited to, seismic waves. The present system and method is directed to 
such a need. 

SUMMARY OF THE INVENTION 

[00010] Accordingly, the present invention provides a Dynamic InSAR imaging 
method that includes sampling fast-moving surface displacement waves from a 
plurahty of moving platform positions and reconstructing a plurality of phase 
differentials from the plurality of platform positions to produce a series of 
interferometric images of said waves. 

[00011] Another aspect of the present invention provides a Dynamic InSAR 
imaging method that includes sampling fast-moving surface displacement waves 
from a plurahty of moving platform positions and subtracting a phase difference 
from phase centers having a time lag tj^g at each of the platform positions. A plurality 
of phase differentials from the plurality of platform positions thus produces a series 
of interferometric images of the fast-moving surface displacement waves. 

[00012] Another aspect of the present invention provides a Dynamic InSAR 
imaging method that includes sampling fast-moving surface displacement waves 
from a plurality of moving platform positions and subtracting a phase difference 
between corresponding antenna formed from offset phase centers having a time lag 
ti3g at each of the platform positions. A plurality of phase differentials from the 
plurality of platform positions thus produces a series of interferometric images 
having one or more sets of concentric rings associated with the fast-moving surface 
displacement waves. 
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[00013] A further aspect of the present invention provides an InSAR apparatus 
that is arranged as a dual phase-center antenna allowing for a phase center offset for 
imaging fast-moving surface waves. 

[00014] Accordingly, the present invention provides an InSAR method and 
system that provides real time velocity assessment of fast moving surface 
displacements of greater than about lOOm/sec, more often between about 300 and 
about 8,000m/sec. Applications include, but are not limited to, earthquake rupture 
physics, volcano harmonic tremor source characterization, open ocean tsunami 
tracking and early warning, and detection/monitoring of damage assessment of 
surface and imderground activity in military and non-military events, such as, 
explosions (e.g., mining). 

BRIEF DESCRIPTION OF THE DRAWINGS 
[00015] The accompanying drawings, which are incorporated into and form a part 
of the disclosure, illustrate an embodiment of the invention and, together with the 
description, serve to explain the principles of the invention. 

Figure 1(a) illustrates an InSAR Strip-map mode geometry. 

Figure 1(b) illustrates a simplified schematic of two SAR images formed 
from two physically separate antennas or a single antenna having two channel 
feeds. 

Figure 2 shows the general geometry for InSAR. 
Figure 3 illustrates a wave sampling schematic. 
Figure 4 illustrates a single antenna design having offset dual-phase 
centers routed to dual receivers. 
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DETAILED DESCRIPTION OF THE INVENTION 
[00016] Referring now to the following detailed information, and to incorporated 
materials; a detailed description of the invention, including specific embodiments, is 
presented. 

[00017] Unless otherwise indicated, numbers expressing quantities of ingredients, 
constituents, reaction conditions and so forth used in the specification and claims are 
to be understood as being modified by the term "about". Accordingly, imless 
indicated to the contrary, the numerical parameters set forth in the specification and 
attached claims are approximations that may vary depending upon the desired 
properties sought to be obtained by the subject matter presented herein. At the very 
least, and not as an attempt to limit the application of the doctrine of equivalents to 
the scope of the claims, each numerical parameter should at least be construed in 
light of the number of reported significant digits and by applying ordinary rounding 
techniques. Notwithstanding that the numerical ranges and parameters setting forth 
the broad scope of the subject matter presented herein are approximations, the 
numerical values set forth in the specific examples are reported as precisely as 
possible. Any numerical value, however, inherently contains certain errors 
necessarily resulting from the standard deviation found in their respective testing 
measurements. 
General Description 

[00018] Conventional Interferometric Synthetic Aperture Radars (InSARs) are 
well understood and extensively used in static mode for imaging static deformations 
associated with earthquake rupture, fault creep, volcaruc deflation, and other static 
deformations from remote targets, but their extension to real-time dynamic imaging 
mode, to image surface displacements produced from, for example, seismic waves 
from earthquakes, explosions, , volcanic generated harmonic tremors to pinpoint 
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magma and gas movement. Tsunamis, storm surges, and oil slick real-time tracking, 
have not been developed. In addition, imaging man-made explosions can lead to the 
detection and characterization of imderground burikers and cavities, which has 
many military uses. The major technical challenges include the need to observe 
small amplitude surface oscillations (often only about 1 cm) that are moving fast 
(greater than about 100 m per second). To form such images, any InSAR platform, 
such as a satellite, or a jet aircraft, must fly in a straight line or smooth arc. If not, 
deviations from a straight line or arc must be measured precisely and compensated 
for during processing (e.g., motion compensation). 

[00019] The present invention. Dynamic Interferometric S5mthetic Aperture Radar 
(Dynamic InSAR) is analogous to "Along-track'' InSAR (discussed herein before, 
and thus far used exclusively to measure relatively slow-moving water current 
waves), but incorporates a sampling method and an improved single anterma design 
that can form real-time images and measure fast-moving surface waves on water 
(Tsunamis) and land (seismic surface waves). Such a system and method can operate 
in the repeat-pass mode or single-pass mode (also discussed above) and within the 
microwave frequency range from about 1mm to about 1 m. Often the present 
invention operates in predetermined frequency bands, such as C-band, e.g., between 
about 2.5 and about 3.5 cm, L-band, e.g., between about 15 and about 30 cm, and S- 
band e.g., between about 8 and about 15 cm. More often the present invention 
beneficially operates within X- frequency band, for example, between about 2.5 and 
about 3.5 cm. 
Specific Description s 
InSAR Geometry 

[00020] The present invention can operate in a Strip-map mode. In such a mode, 
an antenna pointing direction is kept fixed while image data are collected. A time 
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window over which back-scattered radiation is collected using a series of pulses 
from an illuminated target area determines the length of the strip. The collected 
radiation is processed and an accumulation of data is used to construct a 
synthesized image of a target area. To form images, the platform is required to fly in 
a straight line or arc, or, if not, deviations from a straight line or arc require 
corrections during processing. 

[00021] Fig. 1(a) illustrates such a concept and is designated by the reference 
numeral 100. A single antenna (not shown), capable of being divided into a pair of 
chaimels, or two physically separate antennas 104, and 106, as shown respectively in 
Fig. 1(b) can be mounted, for example, along the front fuselage of a moving airborne 
platform 102, such as an aerial jet aircraft, or it can be installed on a space platform 
such as a satellite, to form an "Along-Track InSAR'', as discussed above. Each single 
antenna (not shown) or if arranged as a pair of anteimas 104, 106, as shown in Fig. 
Kb), is capable of measuring a series of return phases for each illumination pulse 
(illustrated as an illumination source beam 108). A beam footprint 110 
approximately equal to XR/d (with X as the illumination wavelength, R is the range 
to target, and d is the length or diameter of the antenna), is created by beam 108 and 
carried along a strip 112 having a length 114. Each pulse is emitted as a chirp 
function (a pulse with a linearly-varying frequency), and reflections from the pulse 
are received back at the antenna and sampled and combined with a replica of the 
outgoing pulse function. Each pulse thus forms one line of along-track SAR image 
data. When multiple pulses are received one after the other in the along-track 
direction, a 2-dimensional SAR image can be formed. Two such SAR images thus 
formed from similar viewing geometries can be combined interferometrically to 
create an interferometric phase map of the ground region. Instruments such as 
gyroscopes can be used to keep the pointing of illununation beam 108 stable. 
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[00022] In another mode of operation, the present invention is capable of being 
operated in Spotlight InSAR, which includes imaging a target by directing an 
antenna at the center of the imaging scene during the entire pass of a space or an 
airborne platform to provide a longer illumination time, (i.e., a spotlight aperture). 
[00023] Such a similar architecture and method is described in a textbook by Fitch, 
J. P., Synthetic Aperture Radar, Springer-Verlag New York Inc., 1988. p. 125-128. TK 
6592.S95F58. 

[00024] Fig. 2 illustrates the general InSAR geometry as designated by the 
reference numeral 200 for cross-track, static InSAR. Two separate passes of the same 
anterma, or two different anteimas (from two different satellites) are denoted as 
imaging positions A^202, and A2204. The distance between imaging positions is 
called the interferometric baseline and is denoted as B. For Dynamic InSAR, the 
imaging positions are offset along the flight path (into the page — ^not shown) but is 
still referred to as the baseline and denoted as the letter B. Similarly a single 
antenna (not shown) with dual phase centers offset along the flight track (azimuth), 
forms the "Along-Track Interferometer'' with phase center offset in the flight 
direction and referred to as the baseline B. In any case, a pair of slant ranges (p) 216 
and (p + Ap) 220 to the same surface element at (z) is measured. A look angle 0 (the 
angle with respect to vertical downward) and a flying height 224 (h) geometry are 
fixed with a and a' being the angle of baselines B and B' with respect to a horizontal 
222. Dashed lines show paths for a second interferogram pair (of the same site 
formed by antenna A^ 202 and A/ 206) for collecting unperturbed topographic phase 
to be removed from the dynamic (seismic wave) collection pair. The path of a 
moving aircraft or satellite is into the page. 

[00025] The measured phase at each point of the two radar images is the sum of a 
propagation part (proportional to the round-trip distance to the on the ground ) and 
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a scattering part due to the interaction of a radar wave with the ground. If the two 
path lengths are p and 6p, the phase difference (() is: 

1. (^--—Lp. 

or 27i/X times the round trip distance difference. 

Using the law of cosines and approximations, Ap can be shown as: 

2. b.p-B\\, 

With B II = Bsin(0 - a) being the component of the baseline parallel to the look 

direction. Combining equations (1) and (2) results in: 

Alt 

3. <^ = T^"' 

which depends on 0, which in turn depends on topography. 

If a second interferogram (denoted by primes and dashed lines) is 
acquired over the same area so that 9 and p are unchanged, the phases of the two 
interferograms can be compared to obtain: 

4. <Z)/<Z)' = 5||/J5|r. 

The ratio of the phases of the two interferograms is equal to the ratio of the parallel 
components of the baselines, independent of topography. 

If a surface displacement occurs, such as, for example, because of an 
earthquake , a volcanic eruption, or man-made explosions, such that each resolution 
element (pixel) has been displaced between observations, the correlated radar 
echoes of the present invention can observe such displacements. 

In addition to the phase dependence on topography, there is also a phase 
change due to the radar line of sight component of the displacement Ap. If the 
displacement occurs in the primed interferogram, the phase ^' is: 

5. <Z)' = (4;^//l)(B||'+Ay9). 
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If the data from the initial tinprimed interferogram are scaled by the ratio 
of the parallel components of the baseline and subtracted from the primed 
interferogram, then: 

Now II y 5 11= fi 'sinC^ -a') IB sin(^ - a) depends on both the illumination 
geometry as well as the topography at each point in the radar image. To evaluate 
equation (6), the topography 6 is subtracted out by using a digital elevation (DEM) 
map of the image area or alternatively, the topography can be subtracted out by 
using the inteferograms themselves. 

The relative sensitivity of the phase measurement to topography and 
displacement can be found by taking the derivative of equation (5) with respect to z, 
as shown in Fig. 2, with z = h - pcosB, for topography and with respect to Ap for 
displacement. 

Using dz = psinGdG, the derivative of equation (5) with respect to topography is: 

7. df = {AnlX){Bcos{d-a)de\ 
and, 

8. dipldz^ {An I Ji){B cos{0 -a)l p sin(^). 

The derivative of equation (5) with respect to displacement is: 

9. d(^ldl^p^{AnlX\ 

Since the distance to, for example, a satellite to the ground p 216, as 
illustrated in Fig. 2, is much larger (e.g., 10s of kilometers) than the distance between 
antenna positions, which is the baseline B between points 202 and 204, the 
right hand side of equation (8) is much smaller than that of equation (9). Therefore, 
the measured phase is much more sensitive to displacement than to topography. 
Consequently, the radar interferometry is more sensitive to a change in topography 
with time than it is to the topography itself. As an example, for a wavelength of 5.6 
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cm (ERS-1 satellite), 1 meter of topography gives a phase signature of 4.3 degrees, 
while for the same pass pair, a 1 meter surface displacement gives a phse signature 
of 12,800 degrees (approximately 3000 times more sensitive). 

Accordingly, the techiuque derived above can be used to measiare surface 
displacements down to about a millimeter level. 
InSAR sampling strategy of the present invention 

A. (single antenna system/repeat pass interferometry for measuring dynamic 
motions) 

Having one antenna requires sampling of target surface displacements 
(e.g., fast moving surface waves of at least greater than about lOOm/sec, more often 
between about 300 and about 8,000m/ sec and having an amplitude of at least about 
1 cm) that can be produced by, for example, earthquakes, mining and gas line 
explosions, military explosions, etc., as a transient but assumed stationary (i.e., over 
the sampling interval during the wave capture pass) perturbation in topography of 
the imaged region. One pass captures an image of the wave (wave capture pass), 
and another topographical repeat pass with less than about a 1 degree angle for 
interferometry collects the unperturbed topography. 

[00026] Fig. 3 illustrates a single antenna repeat-pass mode sampling strategy of 
the present invention and is generally designated as reference numeral 300. Fig. 3 
illustrates that the most optimal portion of a detected wave 302 is the top /crest V4 
region 304 of wave 302 and the bottom/ trough V4 wave region 306 of wave 302, with 
the combined M wave region 308 and 310 in-between top /crest 304 Va region and 
bottom/ trough 306 V4 wave region averaging to zero. A final reconstructed image 
appears as sub-fringes associated with such crests and troughs of the waves where 
the phase difference between passes is directly proportional to the perturbation in 
the local topography due to surface displacements (e.g., seismic waves) traveling 

-12- 
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through the field of view. Sampling calculations are based on this V4 wavelength 
sampling strategy for imaging surface waves originating from a surface or from an 
underground source, such as an explosion- 
Example: 

Calculated imaging parameters: 

PRF (Pulse Repetition Frequency) = 2 x 10 Hz 

Vpiai (horizontal SAR platform velocity (satelUte)) = 8 x 10^ m/sec 

Vpj^j (horizontal SAR platform velocity (airborne)) = 220 m/sec 

Useis (horizontal seismic surface wave velocity) = 3 x 10^ m/sec 

^seis (frequency of seismic surface wave) = 3 Hz 

Keis (wavelength of seismic surface wave) = 1 km 

^seis (period of seismic surface wave) = 0.333 sec/cycle 

Atpp^p(time between each illumination pulse) = 0.0005 sec/pulse 

Therefore, the sampling (i.e., number of output pulses to image) required 

based on imaging parameters shown above and based on a wavelength sampling 

strategy is as follows: 

Ti/4seis (period of crest or trough) = T^^^ /4= 0.08325 sec 
Ti/4seis/^*pRF (* of pulses/^4 wavelength) = 166.67 pulses per Va wavelength 
Vpiat * AtpRp (platform distance traveled per pulse (sateUite)) = 4m/pulse 
Vpiat * ^*PRF (platform distance traveled per pulse (airborne)) = O.lm/pulse 
Accordingly, The Synthetic Aperture Length for each V4 wavelength 

(SAL^,, J is: 

S ALj^^^^j^ (satellite) = (166.67 pulses per M wavelength) * (4m/pulse) = 
668.68 m per V4 wavelength; and 

SAL^4^^.g (airborne) = (166.67 pulses per V4 wavelength) * (O.lm/pulse) = 
18.33 m per wavelength 
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B. (dual channel, dual phase center ATI '^Along-Track Interferometer'' 

system/single-pass ATI mode) 

[000271 Having two channel feeds to a single anterma with dual phase centers (or 
two physically separate antermas) offset along the flight track (azimuth) direction 
can form an Along-Track Interferometer (ATI) as discussed herein before. An ATI 
exploits the short time delay between two similarly collected SAR images from dual 
phase centers spaced in the along-track (azimuth) direction. The radar backscatter 
received at the two antennas are identical except for a change that occurs in the 
ground scattering characteristics dujring a lag time interval (tj^g), i.e., the time it takes 
for the trailing (aft) phase center to travel to the position of the leading (forward) 
phase center. Therefore, a lag time interval for this embodiment is equal to the ratio 
of a baseline B (the distance between the two phase centers) and a SAR platform 
velocity (Vpj J, i.e., t^^^= ^/^piav Subtracting the phase between the two images 
results in a (short) temporal difference occurring during t^^^ being non-zero. 
[00028] Accordingly, an ATI is beneficial in imaging dynamic waves. The present 
invention uses such a technique that is capable of imaging surface waves traveling 
with a velocity greater than about 50 cm/sec (e.g., water perturbations), however 
the present invention is beneficial in imaging fast moving surface waves of greater 
than about 100 m/sec, more often between about 300 m/sec and about 8000 m/sec 
(i.e., 8 km/sec), utilizing a baseline B, for example, of less than about 100 cm. The 
phase difference ^ between the dual phase centers in the ATI system of the present 
invention is as follows: 

10. (p = {47rB /A)(Useis sin(< inc.)IVpiat\ 

where B is the ATI antenna spacing, Vpj^^ is the platform velocity, X is the radar 
wavelength, and < inc., is the radar incidence angle. Assuming a radar incidence 
angle of about 30 degrees and solving for baseline B: 
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11. !B = (6W / InWplax I Useis), 

An interferometric phase ambigxiity arises if the phase difference between 
individual SAR images exceeds tt, which requires that the phase difference (j) be less 
than ;r radians. Accordingly, a beneficial maximum baseline B, can be calculated and 
implemented into the architecture of the present invention for a given set of imaging 
parameters so as to not introduce such a phase ambiguity. 
Example: 

Calculated imaging parameters: 

PRF (Pulse Repetition Frequency) = 2 x 10^ Hz 
Vpiat (horizontal SAR platform velocity (satellite)) = 8 x 10^ m/sec 
Vpj^j (horizontal SAR platform velocity (airborne)) = 220 m/sec 
^eis (horizontal seismic surface wave velocity) = 3 x 10 m/sec 
^seis (frequency of seismic surface wave) = 3 Hz 
^seis (wavelength of seismic surface wave) = 1 km 
^seis (period of seismic surface wave) = 0.333 sec/cycle 
AtpRp(time between each illumination pulse) = 0.0005 sec/pulse 
<inc. (radar incidence angle) = 30 degrees 

Therefore, a critical (i.e., a maximum) baseline B, for various radar 
example wavelengths is as follows: 

Bx (X-band; = 0.03m) = 4.0 cm (Space) and 1.1 mm (airborne vehicle) 
B^ (C-band; = 0.05m) = 6.7 cm (Space) and 1.8 mm (airborne vehicle) 
(L-band; = 0.24m) = 32.0 cm (Space) and 8.8 mm (airborne vehicle) 
Bp (P-band; Xp = 0.69m) = 92 cm (Space) and 25.3 mm (airborne vehicle) 
[00029] Fig. 4 shows an exemplary dual-phase-center antenna of the present 
invention, designated as reference numeral 400, having columns 402 of antenna 
elements 403. Such an antenna 400 also includes two independent phase centers 406, 
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408, created electronically by alternating on and off outer edge conducting strips 
410, 412, controlled and fed by two separate receiver channels 414, 416. For left 
phase center 406 (controlled by receiver channel 1, 414) to be active, left edge 
conducting strip 410 is on and right edge conducting strip 412 is off. For right phase 
center 408 (controlled by receiver channel 2, 416) to be active, left edge conducting 
strip 410 is off and right edge conducting strip 412 is on. The necessary phase center 
offset is determined from equation (11), as shown above, in addition to the SAR 
center frequency wavelength, the velocities of the antenna platform (aircraft or 
satellite), and the velocity of the fast-moving wave being imaged. Such a phase 
center offset is equal to the required ATI baseline, denoted as the letter B 418 for a 
given wavelength SAR configuration. 

[00030] Accordingly, a series of phase differential maps (interferograrris) by the 
method and system of the present invention is obtained by imaging a plurality of 
intereferometric phase differences along a synthetic aperture in either a repeat-pass 
single-antenna (static) mode or a single-pass (dynamic) mode, i.e., single-antenna 
dual-phase receiver or having dual physically separate antennas). Such a method 
and system is capable of measuring natural (e.g., earthquakes. Tsunamis, etc.) or 
man-made (surface and underground explosions) surface displacements having 
wave velocities greater than about 100 m/sec. For example, as stated herein before, 
earthquake rupture physics, volcano harmonic tremor source characterization, open 
ocean tsunami tracking and early warnings are just some of the applications that the 
present invention is capable of addressing. For all such applications, the precision of 
the individual phase measurements can be improved by increasing the illumination 
microwave power utilized in the present invention and, thereby, the detection 
Signal to Noise Ratio (SNR). 
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[00031] The Dynamic InSAR imaging method described herein can additionally 

assess damage to an underground cavity or structure (target) from an exploding 
ordinance (i.e., battle damage assessment, or BDA). Such a method can image the 
concentric rings associated with the dynamic seismic surface waves (i.e., rings of 
surface displacement or disturbance) as they propagate outward from an explosive, 
earth-penetrater-type ordinance directed at a target (analogous to the concentric 
ripples in a pond propagating outward from a stone or other object thrown into the 
pond), 

[00032] Accordingly, three dynamic InSAR imaging results are possible that 
correspond to three different phenomenological outcomes from the exploding earth- 
penetrator ordinance. The first possible outcome is that the energy from the 
ordinance does not sufficiently mechanically couple into the groimd to generate 
seismic waves outward from the point of impact of the ordinance. In this case, most 
of the energy has simply ''bounced off" the target and can be assumed to not have 
caused damage to the imderground target. 

[00033] In the second case, the earth-penetrator does sufficiently mecharucally 
couple into the ground and produces seismic surface waves propagating outward 
from the point of impact of the ordinance. However, there is no secondary set of 
concentric waves propagating outward from the undergroimd target due to a 
secondary moment release that is generated from a partial or complete collapse of 
the underground target. In this case, the underground cavity or structure is 
assumed to have not been significantly damaged as evidenced by the absence of a 
secondary set of concentric waves (rings) propagating outward from the target that 
can otherwise suggest a partial or complete collapse of the underground cavity or 
structure. 
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[00034] The third case is wher\ two, possibly displaced, sets of concer\tric seismic 
surface waves are observed; one from the exploding earth-penetrator ordinance 
generating seismic waves, and the other from the secondary moment release from 
partial or complete collapse of the imdergroiind structure or cavity generating a 
second, delayed in time, set of concentric seismic waves. Computer simulations of 
seismic source and wave propagation scenarios such as the ones described here have 
been designed and executed by the inventor and show that the secondary seismic 
moment release can be several (between about one and about ten) times the seisnuc 
moment release from the earth-penetrator alone, which illustrates the BDA 
discrimination capability afforded by the application of this technique. 
[00035] The Dynamic InSAR imaging methodology described here can also be 
used to both detect, locate, and determine the gross geometry of an imknown 
underground cavity or structure as follows. If an earth-penetrator-type ordinance is 
detonated near an unknown imderground cavity or structure, and the ordinance 
sufficiently couples mechanically into the groimd to generate seismic surface waves 
(i.e., concentric rings of surface disturbance and or displacement), an imdergroimd 
cavity or structure can generate a visible seismic wave shadow on the leeward (far) 
side of the ordinance impact/ explosion point, and can create diffuse reflections in 
the windward (near) side of the ordinance impact/explosion point. This effect is 
analogous to the effect that a floating stick or log has on ripples propagating 
outward from the impact point of rock thrown into the pond. Computer simulations 
of seisnuc source and wave propagation scenarios such as the one described here 
have been designed and executed by the inventor and show that both leeward 
shadows and windward reflections do in fact occur as a result of a perturbation of 
the earth-penetrator-generated wave field. Furthermore, the gross size and shape of 
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the underground cavity can be determined by the size and shape of the wavefield 
shadow and reflechon created henceforth by the underground cavity or structure. 
[00036] The tactical military deplo)nnent scenario of such Dynamic InSAR system 
applicable to BDA and underground target identification can be a moving aircraft, 
such as, a fighter plane with such an imaging system on board, viewing an attack on 
an enemy target area and providing real-time feedback of the BDA or detection and 
location of an unknown underground target as described above. An exemplary 
application is to augment existing moving target indicator (MTI) SAR system 
displays onboard to have this seismic wave imaging and discriminatory capability. 
For example, instead of simply having a Doppler image showing green dots on a 
cockpit display showing man-made moving targets with their respective velocities, 
such a display can also provide visible green concentric circles indicating the 
information described above for BDA and undergroimd target identification. 
[00037] While the invention may be susceptible to various modifications and 
alternative forms, specific embodiments have been shown by way of example in the 
drawings and have been described in detail herein. However, it should be 
understood that the invention is not intended to be limited to the particular forms 
disclosed. Rather, the invention is to cover all modifications, equivalents, and 
alternatives falling within the spirit and scope of the invention as defined by the 
following appended claims. 
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